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We explore the effect of shape - jagged vs. regular - in the jamming limit of very 
gently shaken packings. Our measure of shape e is the void space occupied by a disordered 
grain; we show that depending on its number-theoretic nature, two generic behaviours are 
obtained. Thus, regularly shaped grains (rational e) have ground states of perfect packing, 
which are irretrievably lost under zero-temperature shaking; the reverse is the case for 
jagged grains (irrational e), where the ground state is only optimally packed, but entirely 
retrievable. At low temperatures, we find intermittency at the surface, which has recently 
been seen experimentally. 

1. Introduction 

Why does shaken gravel pack better than shaken bricks do? We explore this question 
via a one- dimensional model jl!2j : despite its simplicity, it exhibits frustration and slow 
dynamics, features j3] which link the fields of granular compaction [415] and glasses |§J. 
The central issue that we probe is the effect of granular shape; our main finding is that 
irregular /jagged and regular /smooth grain shapes have rather different consequences for 
compaction in the jamming limit. 

That glasses or granular media do not have crystalline ground states is well known, as 
is the fact that their attempts to reach their ground states are 'jammed', i.e., hindered 
by long-range interactions. Our model is based on the following picture of jamming. In 
the absence of holes (grain-sized voids), the only way that granular media can compact 
is by grain reorientation. A disordered grain, in our model, 'wastes space'; that is, it 
occupies a net volume equal to its size, plus that of a partial void PEE]- 1 A reorientation 
of this to an ordered ('space-saving') state 'frees up' the partial void, for use by other 
grains to reorient themselves. This cascade-like picture of compaction has been seen in 
the comparisons with experiment jH] of random graph models [5] of granular compaction. 
Also, as there 0, the response of grains to external dynamics is the local minimisation 
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1 Specifically, in our model, each ordered grain occupies one unit of space, and each disordered grain 
occupies 1 + e units of space, where £ is a measure of the trapped void space for a given granular shape. 
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of void space; the ground states so obtained j 1 112 j resemble much more the random close- 
packed state found in granular systems [TUj than the rather unrealistic crystalline ground 
state obtained in earlier work jjj. 

2. The model: definition and ground states 

Grains are indexed by their depth n, measured from the surface of the column. Each 
grain can be in one of two orientational states - ordered (+) or disordered (— ) - the 'spin' 
variables {a n = ±1} thus uniquely defining a configuration. We posit an ordering field 
h n which constrains the temporal evolution of spin o~ n , such that the excess void space is 
minimised - a constraint which is reasonable in the jamming limit. 

The stochastic dynamics in the presence of a vibration intensity T is defined by the 
transition probabilities: 

w(a n = ± -> a n = =p) = exp(-n/£ dyn =f K/T). (1) 

The dynamical length (or boundary layer |l|2p £d y n is a measure of the extent to which 
free surface effects percolate into the bulk; well beyond this, the dynamics is slow, while 
within it, the free surface still has an effect on the dynamics which are relatively fast. 
The local ordering field h n reads 

h n = em- -m+, (2) 

where m+ and m~ are respectively the numbers of + and — grains above grain n. Equa- 
tion (J2J) shows that a transition from an ordered to a disordered state for grain n is 
hindered by the number of voids that are already above it: in fact h n is a measure of the 
excess void space [TT] in the system. 

In the r — > limit of zero-temperature dynamics jHj, the probabilistic rules (0) become 
deterministic: 

a n = sign h n , (3) 

provided h n ^ (see below). Ground states are the static configurations obeying (jSJ) 
everywhere. A rich ground-state structure is achieved for e > 0, because of frustration j^j, 
whose nature depends on whether e is rational or irrational. We mention for completeness 
that the case e < is a generalisation of earlier work [7j, with a complete absence of 
frustration and a single ground state of ordered grains. 

The rotation number Q = e/(e + 1) fixes the proportions of ordered and disordered 
grains in the ground states: / + = Q, /_ = 1 — Q. For irrational e, (J2J) implies that all the 
local fields h n are non-zero. Qualitatively, irrational values of e denote shape irregularity; 
the above then implies that for such jagged grains, the excess void space is never zero, or 
the packing is never perfect, even in the ground state. It turns out j 1 112 j that the ground 
state is in fact quasiperiodic; the local fields h n lie in a bounded interval —1 < h n < e. 

For rational e — p/q, with p and q mutual primes, Q — p/{p + q), and some of the h n 
can vanish. This means that grain n has a perfectly packed column above it, so that it is 
free to choose its orientation. It turns out j 1 1 12 j that orientational indeterminacy occurs at 
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points of perfect packing such that n is a multiple of the period p + q. 2 Each ground state 
is thus a random sequence of two patterns of length p + q, each containing p ordered and q 
disordered grains. The model therefore has a zero-temperature configurational entropy 
or ground-state entropy £ = ln2/(jo + q) per grain. Qualitatively, rational values of e 
imply a regularity or 'smoothness' of grain shape; one could imagine that regular grains 
would align themselves to fit exactly into available voids where possible, in a ground state 
configuration. This in fact happens, leading to states of perfect packing at various points 
of the column and the observed degeneracy of ground states. 

3. Zero-temperature dynamics: (ir)retrievability of ground states 

Zero-temperature dynamics is a theoretical construct; one assumes that, starting with 
a random array of objects, the limit of zero shaking intensity will cause their ground state 
of packing to be achieved. This is neither obvious, nor, as we will show, correct in general. 

Our analysis shows |1|2| that, under zero-temperature dynamics, jagged grains are able 
to retrieve their unique ground state, whereas for regular grains, the true ground states 
are impossible to retrieve. In the latter case, one finds instead a steady state with non- 
trivial density fluctuations above the ground states, which recall the observed density 
fluctuations above the random close-packed state |12|13j in real granular materials. This 
is intuitively understandable; jagged grains in the limit of extreme compaction can only 
'click' into place in a unique way, while the huge degeneracy of perfectly packed ground 
states for regular grains makes any particular one impossible to retrieve, using a random 
shaking dynamics. 

We recall the rule for zero-temperature dynamics: 

cr n ^ sign h n . (4) 

For irrational e with an initially disordered state, this results in the recovery and ballistic 
propagation of the quasiperiodic ground state, starting from the free surface to a depth 
L{t) ps V(e)t. The velocity V(e) = V(l/e) varies smoothly with e, and diverges as 
V(e) ~ £ for e ^> 1 jll'2j . The rest of the system remains in its disordered initial state. 
When L(t) becomes comparable with £ dyn , the effects of the free surface begin to be 
damped. In particular for t ^> (,d y n/V(e) we recover the logarithmic coarsening law 
L(t) ps £d yn hit, observed in related work |7|9j to model the slow dynamical relaxation of 
vibrated sand |T2*j . 

For rational e, things are more complex; the local field h n in may vanish. We 
choose to update such orientations according to a n — > ±1 with probability 1/2, leading to 
a dynamics which is stochastic even at zero temperature. Here, even the behaviour well 
within the boundary layer £d y n contains many intriguing features, while the dynamics for 
n ^> £d y n is again logarithmically slow jll2j . Focusing on the limit £d yn = °o, our main 
result is that zero-temperature dynamics does not drive the system to any of its degenerate 

2 For e = 1/2, for example, one can visualise that each disordered grain 'carries' a void half its size, so 

that units of perfect packing must be permutations of the triad H , where the two 'half voids from 

each of the — grains are filled by the + grain. The dynamics, which is stepwise compacting, selects only 

two of these patterns, H and — | — . Evidently this is a one-dimensional interpretation of packing, 

so that the serial existence of two half voids and a grain should be interpreted as the insertion of a grain 
into a full void in higher dimensions. 
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ground states. The system instead shows a fast relaxation to a non-trivial steady state, 
independent of initial conditions. In this steady state, the local fields h n have unbounded 
fluctuations as a function of depth |l|2j ; since h n is a measure of excess void space, this 
in turn implies the existence of density fluctuations in our model of shaken sand. 
Fig. ^ shows the variation of these density fluctuations as a function of depth n: 



w: 



(h 2 n )^An 2/3 , A ^0.83. 



(5) 



The fluctuations are approximately Gaussian, with a definite excess at small values: 
\h n \ ~ 1 <C W n . Interestingly, this 'nearly but not quite' Gaussian behaviour of density 
fluctuations in shaken granular media has been seen experimentally jH]; the not-quite- 
Gaussianness was in the experiment interpreted via correlations, which as we will show 
below, are also present in our model. 




Figure 1. Log-log plot of W 2 = {h 2 n ) against depth n, for zero-temperature dynamics with e = 1. 
Full line: numerical data. Dashed line: fit to asymptotic behaviour leading to (J5J) (after jH2j ). 



To reiterate, (J3j) implies that the known ground states of a system of regularly shaped 
grains will never be retrieved even in the limit of zero-temperature dynamics. The steady 
state will, instead, be one of density fluctuations above the ground state. The present 
model, to our knowledge, thus contains the first derivation 3 of a possible source of density 

3 A simple scaling argument explains the observed roughening exponent 2/3. Let h n be the position of a 
random walker at 'time' n. The noise in this fictitious random walk originates in the sites m < n where 



the local field h m vanishes. It is therefore proportional to Y^m=i Prob{/i r . 
condition W% ~ J2m=i 1/W m , yielding the power law J5j|. 



0}, hence the consistency 
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fluctuations in granular media |12|13j . which, here, arise quite naturally from the effects 
of shape. This prediction of complex experimental observations |8|12j is all the more 
startling given its origin from a model of such simplicity. 

We turn now to the issue of correlations. If the grain orientations were statistically 
independent, i.e., uncorrelated, one would have the simple result (h%) = ne, while fl^J) 
implies that (h%) grows much more slowly than n. The orientational displacements of 
each grain are thus fully anticorrelated. Fig. |2] shows that the orientation correlations 
c m ,n = (& m (?n) scale as [U2] 

1 / n — m \ . . 

C m ,n « ?>m,n ~ * [wjvj ' ^ 

where the function F is such that F(x) dx = 1. The fluctuations of the orientational 
displacements are therefore asymptotically totally screened: J2n^m c m,n ~ — c m , m = — 1- 
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Figure 2. Scaling plot of the orientation correlation function c m)Tl for n / m in the zero- 
temperature steady state with e = 1, demonstrating the validity of © and showing a plot of 
(minus) the scaling function F (after jll-j )■ 



Summing up, we see that the (orientational) displacements of regular grains are anti- 
correlated 4 within a dynamical cluster whose size scales as n 2 / 3 . Grains separated by 

4 Correspondingly, from a kinetic viewpoint, these results may be interpreted in terms of the time n 2 / 3 
spent by a walker bouncing back and forth between the walls of a cage, where his steps are consequently 
anticorrelated one with the other. 
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greater than the cluster radius are orientationally screened from each other, i.e., the screen- 
ing length also goes as n 2//3 . Consistently, the order parameter Q n = (a n sign/t n ), propor- 
tional to the ratio of the screening length to the total length, goes as n 2 ! 3 /n ~ n -1 / 3 . 5 

Similar anticorrelations 6 in grain displacements have been observed in hard-sphere sim- 
ulations of shaken powders close to jamming J3|; the observation was that grain displace- 
ments along the direction of vibration were strongly anticorrelated, while transverse to it, 
they were uncorrelated. In the absence of free voids, compaction in the simulations hap- 
pened by the phenomenon of bridge collapse] the vibrations coupled to the longitudinal 
displacements of the grains, allowing upper and lower grains in a bridge to collapse onto 
each other, thus minimising the trapped void space between them. This led to the ob- 
served longitudinal anticorrelations; the low intensities of vibration on a jammed granular 
bed did not allow for a coupling with transverse granular displacements. Given the close 
agreement of these earlier results ^3] with our present ones, we can in hindsight justify 
our choice of a columnar model to model granular compaction in the jamming limit. 

Finally, we point out that the landscape of visited configurations in the steady state 
of density fluctuations has a fractal-like structure |H2j . Fig. El shows that, on whichever 
scale we look, some configurations are clearly visited far more often than others. It 
turns out that the most visited configurations are the ground states of the system (empty 
circles). We suggest that this behaviour is generic: i.e., the dynamics of compaction in the 
jammed state leads to a microscopic sampling of configuration space which is highly non- 
uniform, so that the ground states are visited most frequently. This might be expected 
when the system is constrained to compact. However, it should be noted that under 
continuous shaking, the system cannot rest in a ground state even when one is found; 
it evolves constantly, and thereby generates the observed density fluctuations. It should 
also be noted that the entropy reduction AS ~ n 1 ^ 3 resulting from this fine structure 
is subextensive, and therefore negligible with respect to the free entropy Sfl at = nhi2. 
Our model thus provides a natural reconciliation between, on the one hand, the intuitive 
perception that not all configurations can be equally visited during compaction in the 
jamming limit; and, on the other, the flatness hypothesis of Edwards, which states that 
for large enough systems, the entropic landscape of visited configurations is flat [15J. 

4. Low-temperature dynamics 

We now turn to the investigation of the low-temperature dynamics of the model. Our 
main finding is the observation of intermittency in the position of the boundary, or surface, 
layer; this has recently been observed in experiments of vibrated granular beds [To] . 

For rational e, the presence of a finite but low shaking intensity does not change much; 
the zero-temperature dynamics is in any case stochastic, and low-temperature dynamics 
merely increases the effect of noise. However, for irrational e, low-temperature dynamics 
introduces an intermittency in the position of a surface layer, which separates a quasiperi- 

5 By contrast, when e is irrational, earlier orientations influence all successive ones, as the orientation 
correlations c mn do not decay to zero. The order parameter is Q n = 1 identically while, loosely speaking, 
the screening length scales as n. 

6 Treating n as time in a random walk, these anticorrelations also recall the temporal anticorrelations ob- 
served in recent experiments investigating cage properties near the colloidal glass transition (see e.g. |14p. 
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Figure 3. Plot of the normalised probabilities 2 12 p(C) of the configurations of a column of 
12 grains in the zero-temperature steady state with e = 1, against the configurations C in 
lexicographical order. The empty circles mark the 2 6 = 64 ground-state configurations, which 
turn out to be the most probable (after |H2j ). 



odically ordered region near the surface from a steady state of density fluctuations in the 
bulk. 

This happens as follows: when the shaking energy Y is such that it does not distinguish 
between a very small void h n and the strict absence of one, the site n 'looks like' a point 
of perfect packing. The grain at depth n then has the freedom to point the 'wrong' way; 
we call such sites excitations, using the thermal analogy. The probability of observing an 
excitation at site n scales as U(n) « exp(— 2\h n \/Y); the sites n so that \h n \ ~ Y <C 1 will 
be preferred and thus dominate the low-temperature dynamics. These preferred sites are 
such that nQ is closest to an integer, making them look most like points of perfect packing; 
misalignment vis-a-vis Q thus costs the least. The uppermost excitation is propagated 
ballistically (cf. zero-temperature (irrational e) dynamics) until another excitation is 
nucleated above it; its instantaneous position M{t) denotes the layer at which shape effects 
are lost in thermal noise, i.e., it separates an upper region of quasiperiodic ordering from 
a lower region of density fluctuations (J3J). 

Fig. |U shows a typical sawtooth plot of the instantaneous depth Af{t) , for a temper- 
ature T = 0.003. The ordering length, denned as (A/"), is expected to diverge at low 
temperature, as excitations become more and more rare; we find in fact |1|2| a diver- 
gence of the ordering length at low temperature of the form (A/") ~ 1 / (T| In r|). This 
length is a kind of finite-temperature equivalent of the 'zero-temperature' length £d y n, as 
it divides an ordered boundary layer from a lower (bulk) disordered region. Within each 
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of these boundary layers, the relaxation is fast, and based on single-particle relaxation, 
i.e., individual particles attaining their positions of optimal local packing |7|9j . The slow 
dynamics of cooperative relaxation only sets in for lengths beyond these, when the lengths 
over which packing needs to be optimised become non-local. 
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Figure 4. Plot of the instantaneous depth Af(t) of the ordered layer, for e = $ (the golden mean) 
and r = 0.003. Dashed lines: leading nucleation sites given by Fibonacci numbers (bottom to 
top: F n = 89, F l2 = 144, F 13 = 233) (after [IK] I 



5. Discussion 

In this work, we have tried to explain why jagged grains (such as would be found in 
gravel) pack better than smooth ones (such as bricks), when submitted to gentle shaking. 
Our control parameter is e, the void space wasted by the disordered orientation of a grain; 
we have shown that irrational and rational values of this parameter lead to very different 
effects. 

The ground state of jagged grains (irrational e) is unique and only quasiperiodically 
ordered; that for smooth grains (rational e) is highly degenerate and perfectly ordered. 
This is intuitively obvious; the many rough edges of irregular grains need very special 
orientations to click into a perfect packing, whereas rectangular bricks, for example, can 
be arranged in many ways so that they are perfectly packed. 

The very perfection of the ground states for regularly shaped grains makes them impos- 
sible to retrieve stochastically, even in the limit of zero-temperature dynamics; instead, 
the effect of the latter is to give rise to density fluctuations, as predicted by our model, 
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and observed experimentally The 'rough-and-ready' nature of the ground state of 
jagged grains is by contrast quickly (ballistically) retrievable. Clearly, a sharp distinction 
between neighbouring rational and irrational values of e only makes sense for an infinitely 
deep system; for a finite column made of N grains, the distinction is rounded off by finite- 
size effects. In particular, the characteristic features of any 'large' rational e are no longer 
observed when the period p + q becomes larger than N. 

The density fluctuations seen in the case of regular grains have a slightly non-Gaussian 
nature |B] caused by their (anti) correlations; this is reminiscent of dynamical hetero- 
geneities in strongly compacted granular media ^3], as well as temporal anticorrelations 
in cages [Hj. Also, while the macroscopic entropy ^7] of visited configurations in the 
steady state of density fluctuations is consistent with Edwards' 'flatness' hypothesis [To] , 
the microscopic configurational landscape is very rugged, with the most visited configu- 
rations corresponding to the ground states - as might be expected for compaction in the 
jamming limit. Lastly, the low-temperature dynamics for irrational e leads to an inter- 
mittency of the boundary layer separating quasiperiodic order from disordered density 
fluctuations; for shaking at sufficiently low intensities, it should be possible to test our 
detailed predictions |l|2j for intermittency using irregularly shaped grains. 

Remarkably, many of the above features were obtained at a qualitative level in the 
glassy regime of a much simpler model |7j. On the one hand, this allows us to speculate 
that the shape- dependent ageing phenomena seen there could be retrieved here, i.e., that 
conventional ageing phenomena would only be seen for irregular grains (irrational e). 
On the other hand, it is tempting to ask if the directional causality of the dynamical 
interactions present in this model and the earlier one [7|, could be responsible for their 
qualitative similarity, and thus be a necessary ingredient for modelling 'glassiness' ? 
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